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Abstract: Randomly poled nonlinear crystals are shown to be able to 
emit intense ultra-broadband photon-pair fields with properties comparable 
to those coming from chirped periodically-poled crystals. Their intensities 
scale linearly with the number of domains. Also photon pairs extending over 
intervals with durations comparable to one optical cycle can be generated in 
these crystals. 
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1. Introduction 

Spontaneous parametric down-conversion (SPDC) with its production of photon pairs belongs 
to the most fascinating nonlinear optical effects. Two photons comprising a photon pair can be 
entangled in their degrees of freedom as it was first observed by Hong, Ou, and Mandel JT] 
for temporal correlations. Nonlinear bulk crystals served nearly exclusively as sources of these 
photon pairs for many years. However, many highly nonlinear crystals could not be used due 
to impossibility to achieve natural phase-matching conditions. Here, the concept of additional 
periodic modulation of nonlinear susceptibility as introduced by Armstrong i2J has become 
fruitful and resulted in the invention of poling of nonlinear crystals ^ . Highly nonlinear mate- 
rials in which quasi-phase-matching conditions I?) are met can be efficiently used since then. 
They have allowed the construction of bright and versatile sources of photon pairs. 

Moreover, periodical poling has also allowed to tailor the properties of emitted photon pairs 
using nonlinear domains with variable lengths (chirped periodical poling). Domains of different 
lengths in an ordered structure allow an efficient nonlinear interaction in an ultra-wide spectral 
region extending typically over several hundreds of nm ID |6] 111 [D . It has been shown that 
such photon pairs can posses quantum temporal correlations at the timescale of fs and so can 
be extraordinarily useful, e.g., in metrology (quantum optical coherence tomography 0) or 
quantum-information processing ifTOl . 

Alternatively, domains of different lengths can be ordered randomly. The nonlinear interac- 
tion can be efficient even in this case [stochastic quasi-phase-matching] as studies of the process 
of second-harmonic generation indicate iHl fTTl [121 [131 11411 . Here, we show that these structures 
despite their randomness are able to generate spectrally ultra-wide photon pairs at generation 
rates comparable to chirped periodically-poled crystals (CPPC). This shows that contributions 
from the coherent summation of photon-pair fields coming from ordered domains of differ- 
ent lengths in CPPCs are considerably lower than generally accepted. When properly phase 
compensated photon pairs with correlation times at the fs timescale can be emitted from these 
structures. Moreover, these structures have usually smaller requirements with respect to fields' 
polarization properties and orientation of the nonlinear medium |fT5l|. Also fabrication toler- 
ances are less strict in this case H. Practically, technology of poling produces structures with 
relatively large declinations from an ideal geometry which results in reduction of photon-pair 
generation rates. There is the hope that the fabricated randomly poled structures can even out- 
perform their chirped periodically-poled counterparts. The extension to 2D or 3D geometries 
offers additional possibilities including naturally poled materials llTSl . 

2. Spontaneous parametric down-conversion in poled nonlinear crystals 

Quantum state | describing a photon pair generated in the process of SPDC in a poled crystal 
can be written as follows HI: 



where a two-photon spectral amplitude <t> gives the probability amplitude of emitting a signal 
photon at frequency COs and its idler twin at frequency (B,. It can be obtained in the form: 



assuming cw pumping with amplitude ^p, frequency (Op and wave-vector kp. Creation 
operator al{(Oa) in Eq. ([TJ generates a photon with wave-vector ka and frequency (Oa 




(1) 



<i>{(o,,cOi) = g{co,, cOi)^pF{Ak{co,,cOi))d{co° - o), - «,) 



(2) 



(a = s,i) into the vacuum state |vac). Coupling constant g is given as g{(0s,C0i) = 
X^^^ (0) ^COi-CO,-/ [icTt ^ «.((«()«,■(«,■)] where X^'^^ denotes second-order susceptibility, «a(a)„) 
means index of refraction, and c is speed of light in vacuum. Symbol ISk stands for nonlin- 
ear phase mismatch (A^ = kp — — ki). The stochastic phase-matching function F introduced 
in Eq. (|2]i takes the following form: 

F(A/t) = ^(-1)"-'/ dzexp{iAkz). (3) 

„=1 Jz,,-! 

In Eq. symbol Nl gives the number of domains and n-th domain extends from z — Zn— i 
to z = Zn- Positions Zn of domain boundaries are random and can be described as z„ = z„_i + 
Iq + 5l„ (n = 1, . . . ,Nl, Zq = 0) using stochastic declinations 5l„. The basic domain length Iq 
is determined such that quasi-phase-matching is reached, i.e. /q = T^l^k^, AA:o = Afc(a)[',a)?) 
and co^j stands for the central frequency of field a. Independent random declinations 51,, are 
assumed to obey the joint Gaussian probability distribution P: 

^(5L)=(^^exp(-5L^B5L). (4) 

Co variance matrix B is diagonal and its nonzero elements are equal to \/a^. Stochastic vector 
5L is composed of declinations 5l„\ symbol ^ stands for transposition. 

Integration of the expression in Eq. (O leaves us with the following simple formula (valid 
forA^L> 1): 

Fm^-^,t^{-iy^Mi^kzj). (5) 

This formula can be interpreted such that SPDC occurs only in domains with positive suscepti- 
bility with doubled amplitudes and domains with negative susceptibility j*^^' play only the 
role of a 'linear' filler This elucidates why two neighbouring domains form basic elementary 
units for the interpretation of properties of photon-pairs Q . 

The emitted photon pairs can be characterized by mean spectral density n{cOs,COi) of 
the number of photon pairs. The density n is defined along the formula «(fi)v,a),) = 

((V/|fll(to.,)a,(a).v)fl- (co,-)fl,-(a);)|v/))av where the symbol () av means stochastic averaging over 
an ensemble of random realizations of the crystal. In practice mean values can correspond to 
averaging over states of photon pairs emitted at different positions in the transverse plane of 
a naturally -poled material ifTSl . Or they can just serve as an indicator of expected values for 
individual realizations. The spectral density n in quantum state | reads: 

n{co„(o,) = M^ll^l^^(|77(A^(co,,fi),))P)av5(a)«-a),-ft),-). (6) 

In the considered random structures, the averaged squared modulus of phase-matching func- 
tion F, defined as / d5L,P{5L,)\F {Ak)\^ , is obtained as follows: 



H{5k)[l-H{5kf^+^ 
[l-H{5k)]^ 



; (7) 



5k{(0s,(0i) = Ak{cOs,COi) — AkQ. Symbol c.c. replaces the complex-conjugated term and 
H{5k) — exp[/5A:Zo] exp[— (c7AA:)^/4]. A detailed analysis of the formula in Eq. ^ shows that 
{\F{Ak)\^)av is peaked around 5k = 0. The larger the deviation a the broader the peak. 



126 q 




Fig. 1. Photon-pair generation rate A' as a function of the number Ni of domains for an en- 
semble of random crystals with standard deviation a equal to m (solid line), 0.5 x 10^^ m 
(solid line with A), and 2 x 10"^ m (sohd hne with o)\N = J dcOs J dmin{cOs,a)i). 

We consider CPPC for comparison. In this case, the ordered positions of boundaries are 
given as Zn ~ nk) + C'(" —Nl/2)^Iq, ^' = C/Ako, and ^ stands for a chirping parameter. The 
phase-matching function F'^^^'V can then be derived in the form ||6l: 

- (^) (-^) |"*<^<«'/^'' - ""'f^-"^™' ' 

(8) 

fix) = (\/^/2)(yC^jcZo + dk/y/^) and erf means the error function. It holds that the 
larger the chirping parameter ^ the broader the phase-matching function F'^^^'^{Ak). 

3. Photon-pair generation rates and intensity spectra 

In our investigation, we use spectrally degenerate collinear down-conversion from a poled 
LiNbOs crystal pumped at the wavelength = 775 nm. The signal and idler photons oc- 
cur at the fiber-optics communication wavelength = A,- = 1.55 jim. The crystal optical axis 
is perpendicular to the fields' propagation direction and is parallel to the vertical direction. All 

(2) 

fields are vertically polarized and so the largest nonlinear element x^j, is exploited. The basic 
domain length Iq equals to 9.515 jim. A structure composed of Nl = 2000 layers is roughly 
19 mm long and typically delivers 2x10^ photon pairs per second per 100 mW of pumping for 
ordered positions of boundaries. 

The most striking feature of random structures is that the photon-pair generation rate 
increases linearly with the number Nl of domains (see Fig. [T]i- On the other hand, standard 
deviation a plays the central role in the determination of spectral widths ASs and ASi of the 
signal and idler fields. The larger the deviation a the broader the signal- and idler-field spectra 
Ss and Si (see Fig. |2^ and Fig. I2J3 for C7 — 1^ transformation). This can be understood from 
the behavior of spatial spectrum of the X^'^\z) modulation participating in phase matching 
conditions: the larger the deviation a the broader the spatial spectrum of X^'^\^)- However, the 
spectral broadening is at the expense of photon-pair generation rates (see Fig.|2};). 

Similar behavior as observed in random structures has been found in CPPCs consider- 
ing the chirping parameter ^ instead of the deviation a. Our investigations have revealed that 
this similarity is both qualitative and quantitative (see Fig.|2]l. For any value of the chirping pa- 
rameter ^ there exists a value of the standard deviation a such that the spectral widths ASs and 
ASi equal. Moreover, also the photon-pair generation rates are comparable. This behavior is 
illustrated in Fig. |2]for structures with Nl = 2000 domains. The signal-field spectra Ss in CP- 
PCs are extraordinarily wide (larger that 1 /im) for sufficiendy large values of parameter 1^ (see 
Fig.|2t). The signal-field spectra Ss coming from random structures can have the same widths 
assuming sufficiently strong randomness: the transformation curve between the standard devi- 
ation a and chirping parameter ^ conditioned by equal spectral widths ASs is plotted in Fig.|2}3. 
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Fig. 2. a) Signal-field spectral width AS, (FWHM) as a function of chirping parameter i^, b) 
transformation curve between the standard deviation a and chirping parameter assuming 
equal spectral widths AS,,, and c) photon-pair generation rate A' for chirped (solid curve 
with o) and random (solid curve) crystals and their ratio rt^ (r^ = A'/A''^'^''^, dashed curve) 
as functions of chirping parameter [,\Ni = 2000. 




The photon-pair generation rates for random and chirped crystals are compared in Fig. |2};. 
Photon-pair generation rates A^ in random crystals reach at least 80 % of those of CPPCs. 

Comparison of the signal-field spectra for one typical realization of the random structure 
and a CPPC (see Fig. O reveals that spectra of random structures cover a larger range of fre- 
quencies and are composed of many local peaks. Nevertheless, the averaged spectrum of the 
random crystal has the same FWHM as the considered CPPC (see Fig.|3}'). 
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Fig. 3. Signal-field spectrum 5., for a) one realization of the random crystal and b) CPPC 
(solid curve with A) and mean value for an ensemble of random crystals (solid curve with 
o). Spectra Ss are given as Ss = hcOs J d(0in{c0s,C0j) and are normalized such that one photon 
is emitted; a = 2.3 x 10"^ m, ^ = 1 x 10*" m"^, Nl = 2000. 



4. Temporal correlations 

These ultra-wide spectra allow to generate photon pairs with extremely short temporal corre- 
lations that can be measured in a Hong-Ou-Mandel interferometer Ul. The coincidence-count 
rate R„ in this interferometer as a function of relative time delay t between two photons forms 
a typical dip (see Fig. HI that can be described by the formula: 



R„{x) = 1 



exp{i(OpT) J d(0sexp{-2iC0sT:){\F {Ak{(Os, ct)° - «.()) | 
Ro = Jd(0,{\F {Ak{(o,,(0^p - 0),)) l^v. 



(9) 



(10) 



Entanglement times derived from the width of the coincidence-count dip can be as short as 
several fs for both random and chirped structures. Each realization of the random crystal leads 
to a sharp peak as evident from the profile of the averaged coincidence-count rate Rn in Fig.|4^. 
They mainly differ in amplitudes of oscillations occurring at the shoulders. 




Fig. 4. a) Coincidence-count rate R,, and b) sum-frequency field intensity P"™ as they 
depend on relative time delay T for one realization of the random crystal (solid curve), 
chirped crystal (solid curve with A) and an ensemble of random crystals (solid curve 
with o). In b), ideal phase compensation is assumed and curves are normalized such that 
/"„rfT/^"™(T) = 1; values of parameters are the same as in Fig. [5] 



The behavior observed in the interferometer shows the potential to generate photon pairs 
with wave-packets extending over the period of only several fs provided that a good phase 
compensation is reached. The temporal profile of a wave-packet is given by the fourth-order 
stationary correlation function /<i>(t) defined as 



(vac|£,W(0£^'(f-T)|v/> 



T 



2 



(11) 



where T denotes detection duration. Function gives the probability of detecting an idler 
photon in the instant that precedes the instant of signal-photon detection by T. In experiment, 
function /(I)(t) is obtained from the measurement of intensity P™^ of the sum-frequency field 
generated by two photons mutually delayed by T [/<i>(t) °^ /™™(t)]. Intensities /™™ created by 
photon pairs coming from the considered ideally phase-compensated structures are plotted in 
Fig. and demonstrate the principal ability to detect both photons in a window « 5 fs wide. 
Restricting ourselves to quadratic phase compensation |l6l[l6l[l7l, this window is roughly two 
times wider Quadratic compensation is more powerful in CPPCs, however, the difference is 
not large. This emphasizes the potential of the studied photon pairs as a tool of diagnostics of 
ultra-fast processes in physics, biology or chemistry. 



5. Conclusion 

Quantitative similarity in properties of photon pairs generated from randomly poled and chirped 
periodically-poled crystals has been found. Namely ultra-wide signal and idler fields can be 
emitted from randomly poled crystals. Two photons in a pair can in principle occur together 
in a temporal window that characterizes one optical cycle. Importantly, photon-pair genera- 
tion rates are comparable and depend linearly on the number of domains. Contrary to chirped 
periodically-poled crystals the randomly poled crystals are quite tolerant in fabrication. This 
gives a great promise for the use of randomly poled crystals as standard bright ultra- wide spec- 
tral sources of entangled photon pairs. 
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